BELLE 



KEK Preprint 2005-103 
Belle Preprint 2006-7 



Evidence for B rj'n and improved measurements for B rj'K 

J. Schiimann,^^ C. H. Wang,^^ K. Abe/ I. Adachi/ H. Aihara,''^ D. Anipko,^ K. Arinstcin,^ Y. Asano,^'^ 
T. Aushev," A. M. Bakich,^^ V. Balagura," E. Barberio,i* M. Barbero,*^ A. Bay,^^ I. Bedny,^ U. Bitenc,!^ 
I. Bizjak,i2 S. Blytli,2i A. Bondar/ A. Bozek,^^ M. Bracko/'^-i^ T. E. Browder,^ M.-C. Chang,42 p. Chang,23 
Y. Cliao,23 A. Clien,2i W. T. Chen,^! B. G. Cheon,^ Y. Choi,^^ Y. K. Clioi,^^ A. Chuvikov,32 S. Cole,^^ 
J. Dalseno,^^ M. Danilov," M. Dash,^® J. DragicJ A. Drutskoy," S. Eidelman,^ D. Epifanov,^ S. Fratina,^^ 
N. Gabyshev,! T. Gershon/ A. Go,^! G. Gokhroo,^^ B. Golob,^"' a. Gorisek,!^ H. C. YLa.,^'^ J. Haba/ T. Hara,^^ 
N. C. Hastings,43 ^ Hayasaka,!^ H. Hayashii,20 M. HazumiJ L. Hinz,!^ T. Hokuue/^ Y. Hoshi S. Hou,2i 
W.-S. Hou,23 Y. B. Hsiung,23 T. lijima/^ A. Imoto,^^ K. Inami,!^ A. Ishikawa7 R. ItohJ M. Iwasaki,43 
Y. Iwasaki/ J. H. Kang,'*^ P. Kapusta,^" N. Katayama/ H. Kawai,^ T. Kawasaki,^^ H. R. Khan,'''* H. Kichimi/ 

S. K. Kim,36 S. M. Kim,37 K. Kinoshita,^ R. Kulasiri," R. Kumar,30 C. C. Kuo,2i A. Kuzmin,* Y.-J. Kwon,49 
J. S. Lange,^ J. Lee, 36 T. Lesiak,24 A. Limosani/ D. Liventsev," J. MacNaughton,^ G. Majunider,39 F. Mandl,^ 
D. Marlow,32 T. Matsumoto,'*^ A. Matyja,^* W. MitarofF,^ H. Miyakc,^^ H. Miyata,^^ Y. Miyazaki/^ R. Mizuk," 

G. R. Moloney/** T. Mori,'*'' T. Nagamine,''^ I. Nakamura,"^ E. Nakano,^* M. Nakao,'^ Z. Natkaniec,^* S. Nishida,'^ 

O. Nitoh,46 S. Ogawa,'*" T. Ohshima,*^ T. Okabe,!^ S. Okuno,i3 S. L. Olsen,^ H. Ozaki,'^ P. Pakhlov," 

H. Palka,24 C. W. Park,37 H. Park,!^ N. Parslow,38 L. S. Peak,38 R. Pestotnik,!^ L. E. Piilonen,48 M. Rozanska,^^ 
Y. Sakai,'^ T. R. Sarangi,'^ N. Sato,*^ T. Schietinger,!^ Q. Schneider,!^ C. Schwanda,^ R. Seidl,33 K. Senyo,*^ 
M. E. Sevior,i8 M. Shapkin,*" H. Shibuya,*'^ B. Shwartz,* V. Sidorov,* A. Sokolov,*" A. Somov," N. Soni,30 

R. Stamen,'^ S. Stanic,^^ M. Staric,*^ K. Sumisawa,^^ S. Suzuki,34 F. Takasaki,'^ K. Tamai,^ N. Taniura,^^ 
M. Tanaka,'^ G. N. Taylor,*^ Y. Teramoto,^^ X. C. Tian,3i K. Trabelsi,^ T. Tsuboyama,'^ T. Tsukamoto,'^ 
S. Uehara,^ T. Uglov,** K. Ueno,23 S. Uno,"^ P. Urquijo,!^ Y. Usov,* G. Varner,^ K. E. Varvell,38 S. Villa,!^ 
C. C. Wang,23 M.-Z. Wang,23 Y. Watanabe,"'' E. Won,*" Q. L. Xie,^ A. Yamaguchi,42 Y. Yamashita,^^ 
M. Yamauchi,"^ J. Ying,3i Y. Yusa,'*^ C. C. Zhang,^ L. M. Zhang,35 Z. P. Zhang,35 V. Zhihch,* and D. Ziircher^^ 

(The Belle Collaboration) 

^Budker Institute of Nuclear Physics, Novosibirsk 
^ Chiba University, Chiba 
Chonnam National University, Kwangju 
^University of Cincinnati, Cincinnati, Ohio 45221 
University of Frankfurt, Frankfurt 
* University of Hawaii, Honolulu, Hawaii 96822 
''High Energy Accelerator Research Organization (KEK), Tsukuba 
^Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 
"institute of High Energy Physics, Vienna 
^'''institute of High Energy Physics, Protvino 
Institute for Theoretical and Experimental Physics, Moscow 
^''J. Stefan Institute, Ljubljana 
^'^ Kanagawa University, Yokohama 
Korea University, Seoul 
Kyungpook National University, Taegu 
^''Swiss Federal Institute of Technology of Lausanne, EPFL, Lausanne 
^''University of Maribor, Maribor 
^''University of Melbourne, Victoria 

'"Nagoya University, Nagoya 
'''"Nara Women's University, Nara 
''''National Central University, Chung-li 



^'''National United University, Miao Li 
Department of Physics, National Taiwan Univers 
^^H. Niewodniczanski Institute of Nuclear Physics, Krakow 
^'"'Nippon Dental University, Niigata 
^"Niigata University, Niigata 



2 



^'^Nova Gorica Polytechnic, Nova Gorica 
Osaka Gity University, Osaka 



29 



Osaka University, Osaka 



'^'^Panjab University, Ghandigarh 
'^^ Peking University, Beijing 
Princeton University, Princeton, New Jersey 08544 
^^RIKEN BNL Research Genter, Upton, New York 11973 
'^'^Saga University, Saga 
University of Science and Technology of Ghina, Hefei 
'^'^Seoul National University, Seoul 
"''^ Sungkyunkwan University, Suwon 
University of Sydney, Sydney NSW 
'^^Tata Institute of Fundamental Research, Bombay 
^" Toho University, Funahashi 
^'^ Tohoku Gakuin University, Tagajo 
^"^ Tohoku University, Sendai 
'''^'^ Department of Physics, University of Tokyo, Tokyo 
Tokyo Institute of Technology, Tokyo 
Tokyo Metropolitan University, Tokyo 
'^^ Tokyo University of Agriculture and Technology, Tokyo 
^''University of Tsukuba, Tsukuba 
Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061 

Yonsei University, Seoul 
University of Ljubljana, Ljubljana 

We report evidence for exclusive two-body charmless hadronic B meson decays B rj'n, and im- 
proved measurements of B ^ r]' K. The results are obtained from a data sample of 386 x 10^ 
BB pairs collected at the T(4S') resonance, with the Belle detector at the KEKB asymmet- 
ric energy e+e~ collider. We measure B{B+ ^ r)'-K+) = [1.76^0 j2(stat)l:^ :;^(syst)] x 10"^ and 

B{B° rj'n°) = [2.79lJ:°^(stat)lg:25(syst)]xlO-^ We also report the ratio of '^^^^guZt'K") = 
1.17 ± 0.08(stat)±0.03(syst) and direct CP asymmetries for the charged modes. 
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Information on the two-body charmless hadronic B 
meson decays B ^ rj'n is very limited at present Q. 
Measurements of these decay modes can improve the un- 
derstanding of the flavor-singlet penguin amplitude with 
intermediate t, c and u quarks [3- Theoretical predic- 
tions for the branching fractions cover the range (1- 
17)xl0~^ and (l-8)xl0~^ for the charged and neutral 
decays, respectively 0,13 • Recently the charged decay 
was measured by BaBar 0. In contrast, the channel 
B rj'K has been precisely measured 0, 0, 0. In 
the Standard Model (SM) the decay is believed to pro- 
ceed dominantly via gluonic penguin processes 0, and 
has been evaluated with generalized factorization ap- 
proaches 0, 0, ^1 . The measured branching fractions 
are, however, significantly larger than these expectations. 
This has led to speculations that SU(3 )-single t couplings 
unique to the rj' meson or new physics |l2l Il3j contribute 
to the amplitude. More precise measurements are needed 
to constrain the amplitudes and to distinguish between 
theoretical models. Measurements of ratios of branching 
fractions, which reduce the effects of form factor uncer- 
tainties, are especially useful in this regard. 

Additional constraints can be provided by the direct 

CP asymmetry, Acp = B(f H/Hmg^/) ' '^^^^'"'^ f 

final state and / is its CP conjugate. Direct CP violation 



in the rfn^ mode can be large in the SM ^ while 

a non-zero value for Acp m 5+ rj' may require a 
new physics contribution 

In this Letter, we report evidence of _B ^ jyV, im- 
proved measurements oi B ^ rj'K and a direct CP 
violation search in the charged B meson decay modes. 
The results are based on a data sample that contains 
386 X 10® BB pairs, i.e. 35 times larger than our pre- 
vious dataset |2|, collected with the Belle detector at 
the KEKB asymmetric-energy e+e" (3.5 on 8 GeV) 
collider 0. KEKB operates at the T(45) resonance 
(Vi = 10.58 GeV). 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector, a 50- 
layer central drift chamber (CDC), an array of aero- 
gel threshold Cerenkov counters (ACC), a barrel- like ar- 
rangement of time-of-fiight scintillation counters (TOF), 
and an electromagnetic calorimeter comprised of CsI(Tl) 
crystals located inside a superconducting solenoid coil 
that provides a 1.5 T magnetic field. An iron flux-return 
located outside of the coil is instrumented to detect 
mesons and to identify muons. The detector is described 
in detail elsewhere |l5| . Two inner detector configura- 
tions are used. A 2.0 cm beampipe and a 3-layer silicon 
vertex detector is used for the first sample of 152 x 10® 



3 



BB pairs (Set I), while a 1.5 cm beampipe, a 4-layer sili- 
con detector and a small-cell inner drift chamber are used 
to record the remaining 234 x 10^ BB pairs (Set II) flij . 

Charged hadrons are identified by combining informa- 
tion from the CDC (dE/dx), ACC and TOF systems. 
Both kaons and pions arc selected with an efficiency of 
about 86%. Tighter criteria are applied to the pion can- 
didate in B~^ — > rj'-K^, resulting in an efficiency (kaon 
misidentification probability) of about 77% (4%). 

The rj' mesons are reconstructed via two decay chains: 
r]' — > riTT^TT^ (with r] —> 77) and 77' — > p'^'j. We recon- 
struct tt", t], rj' and Kg candidates using the mass 
windows given in Table ^ In addition, we require the 
following. All photons are required to have an energy of 
at least 50 MeV, photons from 77' in rj' p^^ of at least 
100 MeV. The transverse momenta of tt^ for p° candi- 
ates have to be greater than 200 MeV/c. The vertex 
of the TT+TT^ has to be displaced from the inter- 

action point and the Kg momentum direction must be 
consistent with its flight direction. For i?° v'^'^j "^c 
require \Ko\ = em^eIZ] < ^'^^ (O'^) ^' ^ 
(77' P^j), where £^(71.2) is the energy of the two tt*^ 
decay photons. Similarly, we require \hri\ < 0.85. 



TABLE I: Mass windows to reconstruct intermediate states. 



mode 




mass window (MeV/c*^) 


77 




[118,150] 


±2.5(7 






[550,870] 




77 ^ 77 




[500,570] 


+2.5a/-3.3cr 


' + — 
77 —> rjir n 


(in rj'K) 


[945,970] 


±3.4(7 


77' p^ (in 


V'K) 


[935,975] 


±3(7 


77' r]n'^n^ 


(in t/'tv) 


[950,965] 


±2.5(7 


77' ^ p"7 (in 


ti'tt) 


[941,970] 


±2.5(7 


ft) + - 

Kg TT TV 




[485,510] 


±3(7 



B meson candidates are reconstructed combining an 
77' meson with a pion or kaon candidate. Two kinematic 
variables are used to extract the B meson signal: the en- 
ergy difference, AE = Eb — S'bcam! and the beam-energy 



constrained mass, Afbc = V-^boam/'^^" {Pb/c)^, where 
-Bbeam IS the beam energy and Eb and Pb arc the re- 
constructed energy and momentum of the B candidate 
in the T(4S') rest frame. Events satisfying the require- 
ments A'/bc > 5.2 GeV/c2 and \AE\ < 0.25 GeV are 
selected for further analysis. Around 10% of these events 
have multiple B candidates. Among these candidates the 
one with the smallest xlt^ + X^' is selected, where Xvtx 
is a goodness of vertex fit for all charged particles and 
X^, = [(Af(77') - 7Tt^'/(T^']^, where M(?7') is the 77' candi- 
date mass, rrirj' is the nominal mass of the 77' and ajj' = 8 
MeV/c^ is the width of the 77' mass distribution. 

Several event shape variables (defined in the center of 
mass frame) are used to distinguish the spherical BB 
topology from the jet-like qq continuum events. The 
thrust angle 9t is defined as the angle between the 77' mo- 



mentum direction and the thrust axis formed by all par- 
ticles not belonging to the reconstructed B meson. Jet- 
like events tend to peak near | cos 6t\ = 1 , while spheri- 
cal events have a uniform distribution. The requirement 
I cos 6't I < 0.9 is applied prior to all other event topology 
selections. 

Additional continuum suppression is obtained by us- 
ing modified Fox- Wolfram moments ^3 and the angle 
9b between the flight direction of the reconstructed B 
candidate and the beam axis. A Fischer discriminant 
(T) 0] is formed from a linear combination of cos 9t, 
S± and five modified Fox- Wolfram moments. S± is the 
ratio of the scalar sum of the transverse momenta of all 
tracks outside a 45° cone around the 77' direction to the 
scalar sum of their total momenta. These variables are 
then combined to form an event-topology likelihood func- 
tion Cs iJ~-qq)j where s (qq) represents signal (continuum 
background) . We include the quality of the B fiavor tag- 
ging of the accompanying B meson to improve continuum 
rejection. The standard Belle B tagging package 0] is 
used, which gives the B flavor and a tagging quality r 
ranging from zero for no fiavor to unity for unambigu- 
ous fiavor assignment. The data is divided into three r 
regions. Signal-like events are selected by applying likeli- 
hood ratio TZ = Cs/ {Cg + Cqq) requirements optimized on 
MC events in the three r regions separately. For channels 
with an 77' —f ^"7 decay an additional variable 9h, which 
is the angle between the 77' momentum and the direction 
of one of the decay pions in the p'^ rest frame, is included 
for better signal-background separation. 

The branching fractions are extracted using extended 
unbinned maximum-likelihood fits to two-dimensional 
{AE,Mhc) distributions for the 77' rjir^Ti^ and 77' 
p°7 sub-decays simultaneously. The extended likelihood 
function used is: 



L{Ns.Nb^ 



m 



n 



NsPs{AE,,Ah,^) 



-Y,Nb,Pb,{AE,,M^,^) 



(1) 



where Ns (N^Bj) is the number of signal events (back- 
ground events of source .7) with probability density func- 
tions (PDFs) Ps {Nbj) and the index i runs over the 
total number of events TV = Ns + J2j ^Bj ■ 

The reconstruction efficiencies are determined from 
signal MC samples, using the EvtGcn package with 
final state radiation simulated by the PHOTOS pack- 
age I^J]. The efficiencies are calculated separately for 
both Set I and Set II. The absolute efficiency for Set II is 
typically about 0.5% larger than for Set I (for efficiencies 
averaged over the two sets see Tabled). The signal yield 



is expressed as Ns ~ eiN 



(BB)^ 



B 



^(BB)2 ^' where 



B is the signal branching fraction, and and N^gj^y 
arc the efficiency and the number of BB pairs for Set I 
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and Set II. The numbers of B^B~ and B^B^ pairs are 
assumed to be equal. Correction factors due to differ- 
ences between data and MC are included for the charged 
track identification and photon, 7r° and rj reconstruction, 
resulting in an overall correction factor of « 0.9 

The PDF shapes for each contribution are determined 
by MC studies. The signal shapes for Ai? and Mbc are 
assumed to be independent. We model the signal using 
a Gaussian with an exponential tail (Crystal Ball func- 
tion) plus a Gaussian for AE and a Gaussian with 
an exponential tail for Mbc. 

We consider four types of backgrounds separately in 
the fit: continuum, b ^ c and two types of charmless 
decays. Continuum background is modeled by a first or 
second order polynomial for AE and an ARGUS func- 
tion for Mbc. Charmless B decays and b c back- 
grounds are modeled with smoothed two-dimensional his- 
tograms. The contributions from charmless B decays are 
split into two components, one for the decay with the 
largest contribution and one for all other charmless de- 
cays. For r{ , the dominant mode, which is 
modeled separately is i? ^ for rj'K^ it 
is B ^ p"K°; for B° -> 77^° it is B ^ pp; and for 
B+ ri'n+ it is the B+ 7]'K+ feed-down. The feed- 
down in 5+ VjTx^ is modeled with the same PDFs as 
used for the r/'_R'+ signal, shifted and with a corrected 
width in AE. 

The continuum shape parameters are allowed to float 
in all modes. The signal mean and width parameters are 
free for the kaonic modes. For the i?+ ifir^ mode 
these parameters are fixed to the values obtained from 
the charged kaon mode. For 77' 7777"^ tt" modes our 
background MC studies show that no contributions from 
b c decays arc expected. The sizes of background con- 
tributions other than continuum are constrained to the 
values expected from MC studies. The B^ rj'K^ com- 
ponent in the — > r]'TT~^ decay is fixed by the branching 
fraction of 5+ -q'K^ as measured here and the proba- 
bility of kaons faking pions. A simultaneous fit with the 
branching fraction and the charge asymmetry as com- 
mon parameters is used. The resulting projection plots 
are shown in Fig. ^ The reconstruction efficiencies and 
fit results are given in Table ITU 

We find first evidence for the neutral decay: 

B{B^ ^ r;V°) = [2.79tJ:0^(stat)l[;-J5(syst)] x IQ-^, 

and evidence for the charged decay: 

B{B^ ^ r/7r+) = [1.76t°;g(stat)+°:}^(syst)] x IQ-^. 

The ratio of the branching fractions for charged and neu- 
tral B ~* rj'K decays is found to be 1.17 ± 0.08 ± 0.03, 
assuming equal production of B^B^ and B'^B^ . The 
charge asymmetries for the B^ ifn^ and _B+ rj' 
decay modes, listed in Table^I show no significant devi- 
ation from zero. 



Systematic errors are estimated with various high 
statistics data samples. The dominant sources are the 
uncertainties of the reconstruction efficiency of charged 
tracks (3-4%), the uncertainties in the reconstruction 
efficiencies for 77 mesons and photons (3-6%) and the 
uncertainty of the PDF shapes and parameters (1-9%). 
Other systematic uncertainties arise from the K feed- 
down in i?+ rj'-K^ (« 2%), the differences between 
data and MC for AE and Mbc in B+ ri'TT+ (w 4%), 
the Ks reconstruction efficiency uncertainty (4%), the 
uncertainty of the sub-decay branching fractions as given 
by the Particle Data Group (PDG) Hj (1.5%), the num- 
ber of BB mesons produced (1%), the efficiency differ- 
ences due to signal simulation by different MC generators 
(1.4%), the uncertainty in the efficiency (1%) and the un- 
certainty from particle identification (0.7%). The errors 
are added in quadrature and we find the systematic errors 
to be ±5.4%, ±7.3%, t?:^% and for B+ if K+ , 

B° Ty'ifO, B+ ?7'7r+ and B°^ 77^0 decays, re- 
spectively. For the charge asymmetry, efficiency based 
systematic errors cancel out. We estimate the possible 
detector bias on Ac p from the charge asymmetry of the 
continuum background in the B^ s- rj' sample which 
is obtained simultaneously from the fit. We assign 0.02 as 
systematic error both for i?+ rj' K'^ and 77'7r+. Other 
contributions from fitting and normalization together re- 
sult in a systematic error of 0.003 for _B+ rj'K^. For 
_B+ rj'iT^ the uncertainties from PDF shapes and feed- 
down contributions add up to ILq''^^- 

The significance of the B^ Vj-n^ yield is 3.2cr, which 
is calculated as tr = ^2 ln(£,„ax/£o), where £max and £0 
denote the maximum likelihood value and the likelihood 
value at zero branching fraction, respectively. The sys- 
tematic error is included in the significance calculation. 
For B'^ rj'n'^ the corresponding significance with sys- 
tematics is 3.1cr. 

In summary, evidence for B — > rj'ir with greater than 
3(7 significance is found and improved measurements for 
the charged and neutral B rj'K decays are reported. 
The measured branching fractions of -B — !■ rj' K decays 
supersede our previous results 0] and are consistent with 
the measurements of CLEO and BaBar |^ . No charge 
asymmetry is observed in the decay modes — > rj' 
and B~^ — > rj'7T~^. 
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the Nil for valuable computing and Super-SINET net- 
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tract No. 2P03B 01324, Poland); MIST (Russia); ARRS 
(Slovenia); SNSF (Switzerland); NSC and MOE (Tai- 
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TABLE II: Signal efficiencies (etot) with sub-decay brandling fractions included and averaged for Set I and Set II for rl — » tjtt'^tv^ 
and ri' ^"7, total signal yields Ns, branching fractions B, signal asymmetries Acp, and significances a. The first errors are 
statistical and the second (if given) are systematic errors. 





B+ ^ 7^'K+ 


^ 77' A'" 


B+ 7?'7r+ 




etot{ri' —> T^TT + TT ) [%] 

etct(r?'^ A) [%] 
Ns 

Acp 
a 


4.31 ±0.03 
2.78 ± 0.04 
1895.7 ±59.5 
69.2 ±2.2 ±3.7 
0.028 ±0.028 ±0.021 
> 10 


1.19 ±0.03 
1.07 ±0.04 
515.3 ±31.7 
58.9li^;^ ±4.3 

> 10 


2.84 ±0.03 
2.89 ± 0.04 
39.0 ± 13.2 

1 7ft+0.67+0.15 
^- ' "-0. 62 -0.14 

0.20115;^^ ± 0.04 
3.2 


1.72 ±0.02 
1.72 ±0.03 
35.8 ± 12.7 

7Q+1.02+0.25 
^- 'f-O.ge -0.34 

3.1 




5.22 5.24 5.26 5,28 5.3 5.22 5.24 5.26 5.28 5.3 5.22 5.24 5.26 5.28 5.3 5.22 5.24 5.26 5.28 5.3 

J\/bc [GeV/c^] A/be [GeV/c^] A/bc [GeV/c^] A/bc [GeV/c^] 

FIG. 1: AE and Mbc distributions for B+ -> r;'A"+ (a, e), B" rj' K° (b, f), B+ ^ 77^+ (c, g) and B° -> r/V (d, h), 
respectively, for the region A/bc > 5.27 GeV/c^and -0.12 GeV < AE < 0.08 GeV for B° ri'-K° and -0.1 GeV < AE < 0.06 
GeV for all others. 
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